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Abstract 

The present work performs full bandstructure calculations to investigate the structural effect and 
the transition mechanisms of the second harmonic generation (SHG) susceptibility of the a — LiIOz 
crystal. The anomalous inconsistency of associated experimental data of the SHG susceptibility 
tensor is elucidated to be dominated by the structural effect especially on the topology of O- 
atoms. On the manipulation of the structural effect, the modification of SHG susceptibility using 
an external pressure is simulated. The calculations of SHG susceptibility tensor are completed at 
finite frequencies and the static limit. The comparison with the experiments is also incorporated. 
On the analysis of the transition mechanisms, the inter-band transition is determined to entirely 
dominate the whole SHG susceptibility at the static limit; however, the effect of the intra-band 
motion is revealed to be as important as that of the inter-band transition at finite frequencies. 
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I. INTRODUCTION 



Recently the research of second harmonic generation (SHG) of semiconductors has re- 
vived, mainly due to the boost on technologies of the associated second-order phenom- 
ena: the electro-optic effect and the optical rectification. The former, widely utilized in 
the electro-optic modulators to tune the intensity of optical signals, plays a crucial role in 
switching and routing optical broadband networks of the up-to-date telecommunications. [l| 
The latter, usually viewed as an inverse operation to the former, is an important mech- 
anism for generating pulsed tera-hertz radiations, whose applications include the cellular 
level imaging, the semiconductor and high-temperature characterization, the chemical and 
biological sensing, etc. [2j In fact, the SHG susceptibility of semiconductors dominates the 
efficiencies of both aforementioned effects. Thus, the investigation of its physical insights 
not only guides the modification of frequency doubling in the laser engineering, but also 
inspires the innovation of technologies of the associated second order phenomena. 

The alpha-phase lithium iodate (a — LHO3) crystal is a famous SHG material for the 
merits: frequency doubling of the widely used radiation of Yd:YAG, a high nonlinear coef- 
ficient, a weak temperature dependence of refractive index, a high damage threshold, etc. 



3] 



bre this crystal was commercialized via the long-term studying in experiments. 
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17l . y, Q] Moreover, very recently it has been 



ized to the nano-structural systems for the application of nonlinear optical waveg- 
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2l| However, the published experimental data of the SHG susceptibility tensor 



of this crystal exhibits an anomalous inconsistency, such as the variation form —12 to 10 
pV/m about the c/31 component at static limit. [221 This is in contrast to the consistency 
in measurements of the linear dielectric constant. [22| So far this anomaly is still an open 
question. To addressing on this feature, the current investigation performs full bandstruc- 
ture calculations basing on the implement of First-principles simulations. In fact, according 
to the best knowledge, the present work is the first quantum-level simulation of the SHG 
susceptibility for this crystal. 

The present study reveals the structural effect as one of the dominating factors to the 
aforementioned inconsistency of experimental data for this system. Since the macroscopic 
polarization constituted by the dipoles of unit cells defines the SHG susceptibility, the study 
of this effect is motivated by the variation of dipoles for the considerable disagreement be- 
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tween the reported structural refinements 23|, |24| and the existing metamorphosis of this 
polar crystal |25|. In fact, the present results highlight this structure effect especially 
on the influence of the topology of O-atoms to the a — L1IO3 crystal. This argument is 
analogous to the recently reported ideology that modifies the nonlinear optic susceptibility 
of polar-ordered systems via the field- manipulating structural twist. In order to attempt 
further to take advantage of the structural effect, the current study simulates the modifi- 
cation of the SHG susceptibility of the a — LHO3 crystal via loading an external pressure. 
On the other hand, the organic polymer-based systems, the mainstream of nonlinear optic 
materials in the next generation, [lj are not only usually in the form of the assembly of 
curled fibers, but also with the high-flexibility property fitting in the pressure manipulation. 
Thus, the results of structural effect emphasized in present article is very heuristic to the 
engineering about the SHG susceptibility of those organic materials. 

The resolution on the mechanisms of SHG susceptibility of the a — LHO3 crystal is also 
dedicated in this work. According to the previous analysis of bandstructure of the a — LHO3 
crystal by the ' k • p' expansion 27] , the importance of the portion pertaining to the intra- 
band transition has predicted, so the adequate theory should distinctively incorporate the 
contribution of this mechanism. Hence, the length-gauge formalism developed by Sipe et 



al. 
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32i | is adopted for the full bandstructure calculation, which is estab- 



lished within the framework of the time-dependent perturbation for the interaction of the 
independent particle with the field of long- wavelength. Starting from the remedy of the 
zero-frequency divergence in the previous works 33j, |34| , Sipe and Ghahramani included the 
sum rule of the periodic parts of the Bloch functions 35] to vanish the coefficient of the 



divergent term with respect to the crystals with filled bands. [28( Due to the incorporation 
of the sum-rule, the effect of intra-band motion was explicitly included and extracted to be 
defined as independent terms from the portion of the inter-band transition. At the same 
quantum level to calculate the SHG susceptibility of crystals, other methods include: the 
application of "2n+l" theorem of the action functional of time-dependent DFT 36j, the 



3ond-charge model 371], and the application of the Berry phase in macroscopic polarization 



38j . They are able to yield accurate results to consist with the experimental measurements; 



however, the physical insights about the transition mechanism provided from those ways are 
very limited. 

The components of the inter-band transition, the intra-band transition, and the mod- 
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ulation of the inter-band part by the intra-band transition are respectively calculated in 
wide frequency spectra. Owing to most of all published data at the static limit, the present 
study especially addresses on the simulations at this point, in which the adopted formalism 



32. 



391 ] is specified to analytically satisfy the Kleinman symmetry 40|, |4jJ. At finite frequen- 
cies, the Kleinman symmetry breakdowns that is not suited for SHG of cw laser radiation. 
Thus, the associated experimental works were seldom performed before. However, the linear 
electro-optic effect and the optical rectification are always operated in the finite frequency 
regime, especially at high frequencies. Then the present results at finite frequencies are good 
reference to the researches of the associated second order phenomena. 

Since the strongly dynamic charge-redistribution stimulated in the high-power lasers had 



been emphasized |42j, thus once the validity of the perturbation scheme was suspected due 
to the considerable disagreement with the experimental measurements for semiconductors, 
i.e. GaAs [3^. Later the argument was postulated to resolve the doubt that the strong 
screening effect of the semiconductors deteriorates the fault of ignoring the non-local effect 
in the local-density approximation (LDA) of the density functional theory (DFT) to influ- 



ence the performance of full bandstructure calculation. 43, Q, l^j] Thus, the modification 



of the scissor operator was proposed in order to improve the accuracy of the calculation. 



30l . l39l . 1461 ] However, such a non-local problem in semiconductors is determined to disap- 
pear in the present system, according to the resulting band-gap values consisting with the 
experimental data. In fact, the present results manifests the perturbation method still to 
be credible for simulating the SHG susceptibility of the a — L1IO3 crystal; however, except 
the above non-local issue, to consider the consistency with the experiments needs to further 
incorporate other subtle factors, i.e. the structural effect proposed in the present investi- 
gation. Thus, the conclusions of this article are very worthy to refer as applying the same 
scheme to other wide-gap oxides-series nonlinear inorganic crystals. 

II. COMPUTATION AND FORMALISM 
A. FIRST-PRINCIPLES CALCULATIONS 

The calculations are performed within the scheme of the density functional theory (DFT) 
with the generalized gradient approximation (GGA), in which the exchange-correlation po- 
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tential functional is parameterized by Perdew et al. 
tion, the projector augmented wave (PAW) method 
linearzed augmented plane wave (FLAPW) method 



. Both ways of all-electron calcula- 



implements of VASP 



51 



48l. |49I| and the modified full-potential 
5ol are respectively applied by the 



52] and WIEN2k packages 



53 



54 



55 



56(. To simulate the 



equilibrium of Hellmann-Feynman forces, the PAW method is applied to determine the en- 
ergetically optimal structure parameters, in which the atomic relaxation is achieved by the 
conjugate- gradient scheme. The FLAPW method is used for the full bandstructure calcula- 
tion, due to its derivative quantity, electric field gradient (EFG) , to be useful on the analysis 
of structure refinements. Therein the core and the valence states are respectively calculated 
relativistically and semi-relativistically. The valence states include the 2s-states for the Li- 
atom; the 2s- and the 2p-states for the O-atom, and the Ad-, the 5s-, and the 5p-states for 
the /-atom. The muffin-tin radii are set to be 1.8, 1.95, 1.5 A for the Li-, the I- and the 
O-atom, respectively. Inside each muffin-tin sphere, the expansion of associated Legendre 
polynomials for spherical harmonics of the wave function is truncated at I = 10. The value 
of parameter RK max was set to be 8.0. To improve the linearization in the semi-core regime, 
extra local orbit are incorporated for the low-lying states: the 4d-states of the /-atoms and 
the 2s-statess of the Li- and the O-atoms. 



B. FORMALISM 

The presently adopted formalism at finite frequencies is composed of three components: 
the inter-band transitions X2inter(~ ^uj; uj, uj), the intra-band transitions X^ntrai - ^uj; u, uj), 
and the modulation of inter-band terms by intra-band terms X2mod(~ uj, uj), as follows 



31] 



xT(-2uj;uj,uj) = xZter(-^^,^) + xZU-^^^) + xZo d (-^^^), (1) 



abc I o \ € 3 f G?k r nm{ r ml r ln} r %fnm fml fin i > , 

X2inter {~2uJ ] UJ , Uj) = ~2 / "T^l^nm^ ~[- r + ~, r + ? r J? 

Tl J 47T d [UJin - UJ m l) [UJ mn - 2uJ) (u) nd - UJ) (UJl n - UJ) 

abc ( 9 \ - e f ^ a r b c \ r fnl fl 

X2intra\ *UJ;UJ,UJ) — 2 7 2^ nm iU) mn r nm \r ml T ln ) 



+ 2 / ^nmi^mrv n mV ml' Ini L 2 I \ 2 1 \ 

Tl J All Ujf n {UJi n - UJ) UJ^UJmi - UJ) 

o -y fnm r nm {A mn r mn } fnm r nm \^ m l r ln}( UJ ml ~ UJl 



^mni^mn ~ 2u) LO^J^mn ~ 2u) 



(3) 
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X2mod(-2w;^,w) 



dk 

2h 2 J 4^ 



J nml " 



fn 



UJ 



( UJ nl r l m { r mn r nl} UJ lm r n l{ r lm r mn}) 



J nm L -± nm \i r , 



.} 



ram n 



mra V^mn 



(4) 



It is worthy to note that the formula of equation has corrected the original wrong typing 
of the equation (B3) printed in the reference 3jJ- The result of the imaginary part is yielded 
by the analytic solution of the pole u+i5 for all 5 — > +0. Subsequently, the result of real part 
is determined by the solved imaginary portion via the Kramers- Kronig relation. Here f nm is 
defined as f n — fm, which f n deotes the Fermi occupation function to be 1 (0) for the occupied 



(empty) state. The bracket {r h mn rni} is defined as {^ m n r n/} 



— ( f /yC I (yi 

2 V mn nl ' 



mn^nl 



!. Since the 



momentum matrix element is much easier obtained from the resulting periodic Bloch 
function than the position matrix element r° TO (k), then the latter is given via the following 



relationship 



29 



30] 



r ram(k) 



Pnm (k) 



(5) 



imu nm (k) 

where a finite r^ m (k) only exists for uj n (k) ^ uj m (k). The matrix element A^ m has the 
definition, A" m = v® n — f mm 5 where v® n is the intra-band velocity matrix element. The 
problem with the unphysical divergence previously mentioned by Rashkeev at al 39j, also 
occurs with the present study in computing the r° when the bands n and m are nearly 



39J, here also set r" m = when Tiuj nm < e 



degenerate. Referring to the previous treatment 
with the small cutoff e of 5 x 10~ 4 Ryd to remove the influence of this problem. 



At the static limit, the formulism 
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391 ] analytically fulfilling the Kleinman symmetry 



is adopted, respectively expressing the inter-band and the intra-band transition as follows 



X r(0;0,0) = XwSter(0;0,0) + x& tt (0;0,0); 

— 1 rami ml In I I r , r , r \ . 

-'nml [yJmJnl ~T UJ n Jl m ~r UJlJ mn ), 



x&«.(0;0,o) = ^J 



^ml^ln^mn 



X2intm(0; 0, 0) = ^ f^E 1 



Ah J 8vr 



UJ 

mn 



r r 



+ (r a r c 



+ r° 

1 V ' r, 



where (r" 



/;/»' iiii>/:e ^nm^mn:c ' rwnte^ "mni 



and rfL„.„ is defined 



as 



32J 



„fe AC _i_ „C Ai 
mn mn ' mn L -^i 



mn — rn/n J 



UJ;, 



ml v ln) 



(6) 
(7) 

(8) 
(9) 



The special-point sampling method of Monkhorst and Pack [57] is applied for the 
integration of all above formula in Brillouin-zone. To consider the remarkable fluctuation 
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within the k-space to the inter-band energy uj mn and the momentum matrix p nm , to achieve 
the convergence of the integration is only by means of a fine k-mesh. Hence, the samplings 
of 12 x 12 x 11 and the 23 x 23 x 21 meshes were tested for the convergence in advance. Both 
give a difference of about 10 percent, comparable to the error range in the experiments, 
for the calculations of the equation (jH]) and even better consistency of the calculations at 
finite frequencies. Thus, the former is used in the present calculations to give creditable 
estimations. 



III. ANALYSIS FOR STRUCTURE REFINEMENTS 

The unit-cell parameters of two experimental results and the energetically optimal 
structure, determined by the equilibrium of Hellmann-Feynman forces, are respectively 
tabulated in the TABLE fl] In order to simulate the pressure effect to the SHG susceptibility 
of the a — LHO3 crystal, the distorted structure with a z-dimension reduced by 1 millionth 
of the original length with respect to the above optimal structure is calculated; this 
equivalently simulates to load an average pressure of ~ 3770 Gpa normal to xy-plane to the 
original system according to the calculation. After the shrinking, its status at the static 
equilibrium is simulated by the atomic relaxation, whose parameters are also exhibited in 
the TABLE HI The structure refinements of simulated cases are illustrated in the FIG. [U 
The associated results on the bandstructure, the value of band gap E g , and the principle 
component of EFG tensor V zz , are illustrated in the FIG. [2] and the TABLE [J, respectively. 

The resulting bandstructures of all calculated structures only differ in the inter-band 
energies, subjected to the relative atomic positions in the unit cell, and agree with each other 
in the curvature, due to the similarity in the lattice constants. Especially, the energetically 
optimal structure and the associated distortion almost have a same bandstructure. Thus, 
according to the tabulated atomic positions with respect to the fixed /-atom, the different 
topology of O-atoms is deduced to dominate the change of bandstructure near the edges 
by comparing with the results of structure (A) and (B), since the previous study [27J has 
revealed the highly ionic feature of the Li-atom to cause little participation in the states 
near the Fermi-level. The maximum and the minimum of the resulting band-gap values are 
3.934 and 3.297 eV. The former consists with the resonant onset indicated in the absorption 
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58l . l59l . l60j and the transmission spectrum [5J ; however, the latter is comparatively closed 



to the recently experimental result (« 3eV) [6jJ]. Because of the resulting values of E g 
consisting with the experimental results, the effect of non-local defect of DFT-LDA can be 
ignored in the present calculations. 

The resulting V zz is given by the term I = 2, and m = in the full-potential expansion 



of the FLAPW calculation. In experiments, it is determined by the resonant 
frequency of nuclear spin quadruple moment for the nuclei with a spin quantum number 
I > 1. |62j The associated experimental value of the V zz for the /(Li)-atom is —52.8 (0.075), 
determined by the cited measurements: the resonant frequency 151.27 x 10 6 (36.4 x 10 3 ) 
Hz, and the asymmetry parameter 0.017 (0) in the transition from the nuclear magnetic 
moment 1/2 to 3/2 63] (1/2 to 3/2 [g5|); the quadruple moment -0.789 barn 64] (0.04 



barn 



66j). In general, the results of the I- and the Li-atom agree with the experimental 



values; it implies the accuracy of the present FLAPW calculations, since the EFG itself can 
be viewed as a criterion on capturing the detail electronic structure of the system. 

According to the similarities of the resulting E g and V zz for all simulated structures 
and the agreement with the experiments, these structures should be equally possible to be 
operated in the lab. Thus, the inconsistency of results defined by them can be considered 
as the nature in the experiments. 



IV. SECOND HARMONIC GENERATION RESULTS 



Static Limit 



The resulting independent nonzero components: d^i and d^; the associated experimental 
results are shown in TABLE [III and IIHI The ^omponent^c^ can be directly measured in the 

IQ, 11, 12,[l3,Q; however, the 
121 ] via another crystal, whose SHG 



crystal itself by the phase-matching method 



7. 



10 



8. 



9 



<i 33 has to be measured by the 'wedge' method 
susceptibility needs to be known. In order to clinically reduce the erroneous judgement from 
any external factor, the comparison with the experiments more emphasis on the component 
of d%\ in the present study. Herein the tabulated portion of the visual-hole (visual-electron) 
process in the inter-band transition is defined by the Eq.([7j), in which the transitions com- 



8 



pletes via two valence (conduction) states. On the other hand, due to the lack of consistent 
experimental data to refer for the a — LHO3 crystal, the test on a well-known system for the 
adopted full bandstructure calculation of the SHG susceptibility becomes necessary, which 
can be viewed as a kind of 'alignment'. Thus, the well-studied GaAs system was selected. 
Basing on the same FLAPW method to compute the bandstructure and the scissors op- 
eration suggested by Nastos et al 46] to deal with the non-local defect of DFT-LDA, the 
resulting value 185 pm/V is rather closed to the experimental result 168 pm/V and the pre- 
viously simulated results collected in the reference [39J]. In general, the resulting magnitude 
of the SHG susceptibility of the a — LHO3 crystal is at the same order of the experimental 
data. 

The importance of the structural effect to the SHG susceptibility of the a — LHO3 crystal 
reflects on the significant variation of the results for different structures on both d$i and 0^33 
components. First of all, basing on the previous study on the bandstructure by the ' k • p' 
expansion the agreement on the curvatures of obtained bandstructures for different 
structures illuminates the consistency of the resulting momentum matrix elements Pm n > so 
the influence of p a mn to the structural effect can be excluded. Thus, according to the anal- 
ysis of bandstructures for all structures in Sec. IHfl this structural effect are deduced to be 
caused by the change on the inter-band energies u mn , which is dominated by the topology 
of O-atoms, previously also mentioned in Sec. IIHI Furthermore, the negative result of the 
structure (A) is in contrast to the positive results of the other structures, can be attributed 
to its unique inverse O-topology among all calculated structures, depicted in the FIG. [Q 
This inverse O-topology feature causes a negative dipole of the 10% anion to compose a neg- 
ative macroscopic polarization of the system. Such a polarization is completely equivalent 
to that caused by the presently obtained negative susceptibility in the same system. Thus, 
the present investigation illuminates the structural effect, especially on the topology of the 
O-atoms, dominating the magnitude and the sign of the SHG susceptibility tabulated in the 
TABLE HTl and HTT1 This conclusion can explain why the inconsistency of the data appears to 
the different specimens in spite of repeating measurements via the same operation conditions 



in the lab. 
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14j 



On the respect of manipulating the structural effect via an external pressure, the magni- 
tude of the modification is about 10 percent to the d$i according to the present simulation. 
Although this scheme seems neither efficient on the low enhancement nor economic on the 
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requested pressure being as high as 3770 Gpa; however, it is still an effective way to modified 
the SHG susceptibility of the a — L1IO3 crystal. 

On the resolution of the transition mechanisms, the contribution of the inter-band tran- 
sition distinctively prevails that of the intra-band transition to both components. Moreover, 
the analysis on the compositions of the inter-band transition reveals the weights of the visual- 
hole and the visual-electron process differing between d$\ and c/33 components, even though 
the results of both components given by the same bandstructure. This feature highlights the 
extreme anisotropic character of the SHG susceptibility of the a — LHO3 crystal. In addi- 
tion, in fact the structural effect significantly changes the proportions of the inter-band and 
the intra-band transitions, and those of the visual-hole and the visual-electron components 
in the inter-band transition. 



B. Finite Frequency Regime 

The results of components d^i, 0^33, du, and di$ at the infinitesimal frequency are tabu- 
lated in TABLE HV] \V\ Wl\ and IVIII Basing on the comparison with the experiments to the 
components du and di§, the resulting magnitude in structure (A) and (B) are comparable to 
the measured values. Thus, the present study manifests the full bandstructure calculation is 
valid to yield the results consisting with the experimental data at the infinitesimal frequency 
for the a — LHO3 crystal. In spite this point infinitively neighboring to the static limit, the 
contributions associated with the effect of intra-band motion suddenly raise up in contrast 
to the status at the static limit, reflecting on the differences between the tabulated results of 
dsi in TABLE [TTJ and IIVI The cancellation because of the resulting signs of the intra-band 
and inter-band portions being appositive makes the magnitude of total result at this point 
much less than that at the static limit. The same trends are also happen to the component 
g?33. On the other hand, the aforementioned structural effect is still very remarkable that 
changes the importance of individual transition mechanism in different structures. To the 
analogous results in the calculations of GaAs, the result 115.8 pm/V is also near to the 
corresponding experimental value, 99.8 pm/V (the value of linear electro-optic coefficient), 



according to the relationship 30] x xyz {~ u 'i u i 0) = X xyz {~ 2^; to) for uj to be infinitesimal. 

The resulting different contributions to the imaginary part of the component 
X2 U {~ 2cu; lu, lu), shown in FIG. El elucidates the intra-band contribution to be as import 
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as the inter-band one in the dispersion at finite frequencies. The structural effect is still 
very distinguished according to the resulting dispersions of different structures. In addition, 
the modification of the pressure effect is appreciable within the low frequency regime but be- 
comes considerable at some high frequency channels to the individual transition mechanism. 
Besides, the previously predicted double-resonance feature as u near to the E g according 
to the analysis of the bandstructure, [27] is reified in the maximal resonant absorption of 
both inter-band and intra-band transitions exhibited in FIG. |3j Especially, this feature very 
insists regardless of the structural effect. On the other hand, the analogous results in the 
calculation of GaAs, exhibited in FIG. El very agree with the published results of the refer- 



ences 



30 



39]. 



The dispersions of all nonzero independent components of the SHG susceptibility tensor: 
xl 23 (—2u] lu, lu), X2 31 (— 2c<j; to, lu), x| n (— 2u;; lu, lu), and X2 33 (~ 2u;; to, lu), are calculated for the 
energetically optimal structure, depicted in FIG. HI The remarkable anisotropic property of 
the SHG susceptibility of the a — L1IO3 crystal reflects on the significant difference among 
the results of those components. 



V. SUMMARY 



The full bandstructure calculation is valid to generate correct results of SHG susceptibility 
to analyze the trends of experimental data at the static limit and at finite frequencies for the 
a — LHO3 crystal. According to the agreements with the experiments of the resulting EFG 
and Eg, the bandstructure give by the FLAPW calculations is valid, so the non-local defect of 
DFT-LDA in the semi-conductors can be ignored in the present calculations. The addressed 
structural effect indeed significantly influences the SHG susceptibility of the a — LHO3 
crystal. Thus, it is proposed as one of the major factors causing the anomalous inconsistency 
of the experimental data collected from published references. On the manipulation of the 
structural effect by an external pressure, the present simulation manifests this scheme to be 
effective. The present work resolves the mechanisms of SHG susceptibility of the a — LHO3 
crystal, and makes the conclusions: the inter-band transition is determined to dominate the 
whole SHG susceptibility at the static limit; the effect of the intra-band motion suddenly 
raise up to be as important as the inter-band transition at finite frequencies. The results 
on the this respect provide very useful information for the further modification of frequency 
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doubling and associated second order phenomena technologies to the a — LHO3 crystal. 
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TABLE I: The parameters of the unit cell, the resulting V zz (in the unit of 10 21 V/m 2 ), and 
the obtained band gap E g (eV) of the simulated structures. The coordinates of each atomic 
species are respectively described as Li-atom(0, 0, z), I-atom(l/3, 2/3, 0), and 0-atom(x, y, z). The 
designations of the structures are as followings, (A): the result of experiment [241 ] ; (B): the result 



of experiment 
(C). 



23J ; (C): the energetically optimal structure, and (D): the pertaining distortion to 



Structure 


o(A) 


c(A) 


Li(z) 


0(x,y,z) 


yLi 
v zz 


V zz E 9 


(A) 


5.48150 


5.170900 


0.8907 


(0.0936, 0.3440, 0.1698) 


0.036 


-53.244 3.735 


(B) 


5.48169 


5.172370 


0.0713 


(0.2468, 0.3419, 0.8377) 


0.057 


-54.016 3.934 


(C) 


5.36816 


5.019063 


0.0775 


(0.2655, 0.3387, 0.8327) 


0.044 


-52.514 3.297 


(D) 


5.36816 


5.019058 


0.0774 


(0.2656, 0.3391, 0.8328) 


0.044 


-52.584 3.304 



[66] H.-G. Voelk and D. Fick, Nuclear Phys., A530, 475 (1991). 
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TABLE II: The d^i component of SHG susceptibility tensor (in the unit of pm/V) at the static 
limit. The calculated results of the structures (A), (B), (C) and (D) and the experimental data 
are tabulated. The 'TOT' labels the result of equation the 'Inter' ('Intra') labels the result 
of inter-band (intra-band) transition from the equation (|7|) (equation (HI)), and the 'V.H.'('V.E.') 
labels the result of the visual-hole (visual-electron) process in the inter-band transition from the 
equation ([7]). The A labels the operated wavelength of the cited experimental works. 



Structure 


Intra 


Inter 


V.H. 


V.E. 


TOT 


(A) 


-0.112 


-3.891 


-4.521 


0.630 


-4.003 


(B) 


0.032 


2.671 


3.330 


-0.659 


2.703 


(C) 


0.031 


4.168 


3.856 


0.312 


4.200 


(D) 


0.014 


4.500 


3.893 


0.607 


4.515 


Experiment 


Result 


A (nm) 








Exp [41 


-12.2 ± 1.9 


1860 








Exp [51 


-4.5 ± 0.6 


1856 








Exp [6] 


-7.54 ± 1.13 










Exp [71 


-4.96 ± 0.26 


1064.2 








Exp [8j 


±10.17 ± 2.0 


514.5 








Exp [9] 


-5.53 


1064.2 








Exp [10] 


-7.215 










Exp [10] 


-7.33 










Exp [10] 


-6.82 


1318 








Exp [11] 


-4.1 


532 








Exp [12] 


-4.4 


1064 








Exp [12] 


-4.0 


488 








Exp [14] 


-4.09 


1064 








Exp [14] 


-3.90 


660 








Exp [14] 


-5.23 


403 









16 



TABLE III: The component of the SHG susceptibility tensor (in the unit of pm/V) at the 
static limit. 



Structure 


Intra 


Inter 




V.H. 


V.E. 


TOT 


(A) 


0.240 


-0.577 




0.623 


-1.200 


-0.337 


(B) 


-0.171 


3.444 




1.724 


1.720 


3.272 


(C) 


-0.104 


1.677 




0.061 


1.616 


1.574 


(D) 


-0.119 


1.801 




0.118 


1.684 


1.683 


Experiment 


Result 


A (nm' 


) 








Exp [5j 


-3.6 ± 1.08 


1719 










Exp [7] 


-5.15 ±0.32 


1064.2 










Exp [10J 


-6.75 ±0.95 


1318 










Exp [10J 


-5.54 ±0.61 


1318 










Exp [12] 


-4.5 


1064 










TABLE IV: The dzi component of the SHG susceptibility tensor (in 


the unit of pm/V) at infinites- 


imal frequency. The 'TOT', the 'ter', 


the 'tra', and the 


'mod' label the contributions of the portion 


defined respectively from the equations (JTJ) , d2J) , 


©, and (H). 






Structure 


ter 




tra 




mod 


TOT 


(A) 


4.084 




-7.812 




3.393 


-0.335 


(B) 


-2.681 




6.409 




-3.917 


0.189 


(C) 


-4.105 




5.152 




-0.168 


0.879 


(D) 


-5.257 




6.325 




-0.189 


0.879 
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TABLE V: The ^33 component of the SHG susceptibility tensor (in the unit of pm/V) at infinites- 
imal frequency. 



Structure ter tra mod TOT 

(A) 0.545 -1.508 -0.712 -1.675 

(B) -3.435 3.204 1.571 1.340 

(C) -1.655 4.800 2.367 5.512 

(D) -1.780 4.817 2.367 5.404 



TABLE VI: The di4 component of the SHG susceptibility tensor (in the unit of pm/V) at infinites- 
imal frequency. 



Structure ter tra mod TOT 

(A) -0.440 -1.759 1.948 -0.251 

(B) 0.000 0.189 -0.398 -0.209 

(C) -0.042 3.540 1.990 5.488 

(D) 0.335 3.016 2.115 5.466 

Exp. [15,16J 0.22-0.35 (A= 1064.2 nm) 



TABLE VII: The c?i5 component of the SHG susceptibility tensor (in the unit of pm/V) at in- 
finitesimal frequency. 



Structure ter tra mod TOT 

(A) 3.854 1.738 -2.618 2.974 

(B) -2.639 -2.807 2.157 -3.289 

(C) -4.147 1.299 2.115 -0.733 

(D) -4.440 1.864 1.843 -0.733 
Exp. [17] 5 (A= 694.3 nm) 
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FIG. 1: (Color online) The different structure refinements of the a — LHO3 crystal. The structure 
(A) and (B) are given by the x-ray diffraction results 



2| and 



231 ] . sketched as left and right figures, 



respectively. The refinements of the energetically optimal structure and the distorted structure in 
simulating the effect of an external pressure along the z-axis, are similar to that of structure 
(B). The dipoles of the IO^ 1 anions are sketched by the green arrow symbols, where the ones of 
structure (B) are defined as positive. The /-atom, the O-atom and the Li-atom are sketched as 
the purple, the red, and the grey spheres. 
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r m k r a 

FIG. 2: (Color online) The bandstructure for the simulated structures. The blue, the red, the 
black, and the green lines respectively depict the results of the structure (A), (B), (C), and (D). 



20 



E 



a 

o 

o 













i — 1 — i — S — i — 4 — t — i — I 










V (A) - 



2 3 4 5 6 7 8 9 10 11 12 13 14 15 

E(eV) 



-1 

















2 3 4 5 6 7 8 9 10 11 12 13 14 15 

E(eV) 



A i 









2 3 4 5 6 7 8 9 10 11 12 13 14 15 
E(eV) 



£V A ft 1 / 






/ (D) ; 



1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 
E(eV) 



FIG. 3: (Color online) The different contributions to the imaginary part of the component 
X2- 1 (—2uj;uj,uj). The results (in the unit of 10 esu) of the structure (A), (B), (C), and (D)are 
exhibited. The results (in the unit of 10 -6 esu) of the GaAs is depicted as the inset of panel (A). 
The contributions of the inter-band contribution, the intra-band contribution, and the modulation 
of the inter-band portion by the intra-band transition, are depicted as the blue, the red, and the 
green lines, respectively; the total is represented by the black line. 



21 




1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 
E(eV) 




1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 



E(eV) 




1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 
E(eV) 

FIG. 4: (Color online) The dispersions of all nonzero independent components of the SHG suscep- 
tibility tensor (in the unit of 10 esu). The black, the red, and the blue lines depict the absolute, 
the real, and the imaginary results, respectively. 
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